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Many-particle entanglement is a striking feature of 
quantum many-body systems. Entanglement was first 
recognized as a curiosity of quantum mechanics because 
it gives rise to seemingly nonlocal correlations of mea- 
surement results of distant observers. Whereas the cen- 
tral role of many-body entanglement for various applica- 
tions of quantum information processing (e.g. (H) is un- 
doubted, its role in e.g. quantum phase transitions (e.g. 
0) or ionization processes is stiU debated (e.g. Q), and 
questions concerning e.g. its potential assistance to the 
astonishing transport efficiency of biological compounds 
(e.g. [4|) are still essentially open. 

To answer such questions we need reliable techniques 
to characterize entanglement properties of general quan- 
tum states. However, even the conceptually rather sim- 
ple question 'Is a given quantum state entangled or not?' 
is in general unanswered so far. It is usually addressed 
by means of separability criteria, which work very well in 
many cases, but are far from perfect Even more chal- 
lenging is the detection of genuine multipartite entangle- 
ment, which has already been intensely studied (see for 
example 0-01 ), but still has not yielded satisfying results. 
Vast areas of the considered state-spaces are still widely 
unexplored due to the lack of suitable tools for detecting 
and characterizing entanglement. 

The central difficulty arises from the complicated 
structure of multipartite entangled states: even states 
that do not separate into blocks of subsystems that are 
not entangled with each other are not necessarily gen- 
uinely n-body entangled. Recently, inequalities to iden- 
tify genuinely n-body entangled states have been pro- 
posed based on nonlinear functions of matrix-elements 
[lOj. Although these new criteria are promising in the 
sense that they allow us to characterize states as entan- 
gled that can not be detected with the standard criteria, 
it is also evident that the characterization of entangled 
states will not be facilitated by a huge set of separabil- 
ity criteria unless we have a systematic way to construct 



and understand these criteria. Here, we present a very 
general, systematic approach to construct such criteria, 
and show that our newly constructed criteria are stronger 
than all formerly known ones. In particular, all these cri- 
teria apply to systems of arbitrarily many subsystems of 
arbitrary finite dimensions. 
In more detail, we derive 

- an TO-linear inequality ([3]) and its bilinear version 
ineq.(jl| to detect bipartite entanglement. Based on this, 
we derive 

- a general framework to obtain bilinear inequalities 

which characterize genuine multipartite entangle- 
ment and 

- construct a particularly strong criterion, i.e. ineq. (|III|) . 
for which the efficiency is demonstrated in the consecu- 
tive examples. 

A pure n-partite state l^*) is called k-separable if it 
can be written as a product [5| 

1^-) = |0l) (g) 102) ■ • ■ » |<?!>fc) , (1) 

of k states each of which corresponds to a single sub- 
system or a group of subsystems. If there is no such form 
with at least two factors, then jvl/) is considered genuinely 
n-partite entangled. On the level of pure states the ques- 
tion of /c-separability can be answered in a straight for- 
ward fashion by means of separability criteria for bipar- 
tite systems, simply by considering all segmentations of 
the fc-partite system into two parts. However, the same 
question becomes significantly more difficult to answer 
for mixed states g: here, a state is considered genuinely 
fc-partite entangled if any decomposition into pure states 

, (2) 

i 

with probabilities pi > contains at least one genuinely 
fc-partite entangled component. Therefore, a mixed state 
can still be partially separable, even if the fc subsystems 
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can not be split into two groups that are not entangled 
with each other. Consider for instance the tri-partite 
state 

Pb^sep = Y1p:'Pab'»Pc+Y1 'ij Pac® Pb+^'' 3 Pbc® Pa ' 

3 3 3 

Here the two-body states p'^^, p'g(j and p'^p describe 
entangled states. Even though there is no bipartite split- 
ting with respect to which the state p is separable, it is 
considered biseparable since it can be prepared through 
a statistical mixture of bipartite entangled states. 

To be certain that some state is really genuinely n- 
body entangled, we thus have to make sure that there is 
no pure state decomposition with only at least partially 
entangled components. Since this reduces to the problem 
of deciding whether each of such pure state components 
is at least biseparable, let us first introduce a suitable 
criterion for biseparability, which then will turn out to 
be the central building block for the subsequent gener- 
alization to genuine many-body entanglement. What we 
employ here, are m-linear functions of a quantum state 
Q on "Ha ®T~iB that can be expressed in terms of the m- 
fold tensor product g®™ of the density matrix g acting 
on the TO- fold tensor product space (%a ^Hb)®™. As it 
is shown at the end of our letter, any separable state 
satisfies 

^|5?e(($|(i ® ^s)t^>f™(^^ ® i)|$))| < ^mgf^\<^) ,(3) 

for any positive integer to, where |$) is any fully sep- 
arable state of the TO-tupled system, i.e. |$) factorizes 
into 2to single-body states. is the cyclic permutation 
operator acting on H^™, i.e. 

Wa\^i)®\(P2)®- ■ ■®|</5m) = |V2)<8)|<^3)'X)- ■ .(E)\(pm)<E)\(pi) , 

(4) 

and Hb is defined analogously for subsystem B. In 
our following extension to multipartite systems, we will 
content ourselves with the bilinear case to = 2, as it 
is already very powerful in detecting entanglement and 
ineq. ([3]) takes the rather simple form 

\{U\p\kj}\-^{zj\p\ij){kl\p\kl)<0 , (I) 

which corresponds to the choice |$) — \ijkl). 
For our following generalization of ineq. ([3]) to the mul- 
tipartite case we will consider all (2"~^ — 1) different 
partitions of an n-partite systems into two subsystems, 
because a mixed state is biseparable exactly if there is a 
decomposition into pure states each of which is separable 
with respect to some partition. The fictitious subsystems 
will be labeled Ai and B,{i = l,..., 2"-^ - 1) in the fol- 
lowing. Introducing the global permutation operator IT 
which performs simultaneous permutations on all subsys- 
tems, we can formulate now the generalization of ineq. ^ 
to multipartite systems: 

^($|p«52n|$) _ \/^MPh^^^Pi\^) < , (11) 



with Vi = IlAi <8 Isi, and where the sum runs over all 
inequivalent bipartitions. 

To convince ourselves that ineq. (|n| is indeed satisfied by 
all at least partially separable states p, let us first verify 
that this holds for any pure state = |4'f,s)(^'fcs| that is 
biseparable with respect some partition labeled io. Just 
like any duplicated state l^')'^^ is invariant under the 
global permutation 11, the duphcated state |^6<,)®^ is 
invariant under II^. (E" 1_b • Therefore, the first term in 
ineq. (|Il|) cancels with the i = io term in the summation. 
All remaining terms are expectation values of positive op- 
erators, and given the negative sign in front of the sum, 
the left-hand-side is indeed nonpositive. Hence, ineq. pi)) 
is satisfied for any pure state that is not genuinely mul- 
tipartite entangled. 

The generalization of ineq. (|ll| to mixed states is a 
direct consequence of its convexity which we can see in 
the following, where we will use that the state 1$) is 
completely separable. That is, independently of which 
decomposition of the Hilbert space into two subspaces we 
take, we can always write it as a direct product of two 
states |$i) and |<1'2) of the respective subspaces. The 
first term in ineq. pll) is the absolute value of the matrix 
element {^i\p\^2)'- 

^(<I,|p82n|$) = |(cOi|p|ci>2)| , (5) 

since 11 simply permutes |$i} and \^2), i-e. n|$i} 
^'2) — 1^2} ^ |$i)- And the absolute value is convex, 
i.e. \a + b\ < \a\ + \b\ for arbitrary complex numbers 

a and b. Each summand K-i = \J {^v} p®'^'Pi\(^) in the 
second term of ineq. (|TT]) is the square root of a product of 
two diagonal density matrix elements, i.e. non- negative 
numbers 

/C, = ^($i|p|$i)($2|p|$2) , (6) 

with |$i) ® |$2) = Ha, ® 1-^*1$). Now, Cauchy- 
Schwarz's inequality J2jPjl3 > \lY.jP]\jY.j 1] with 

P3 = and = y'($2|Pj|$2) yields 

/C. < J2 \/(*l|P.#i>(^2|p,|$2) , (7) 

3 

for any set of positive operators pj satisfying p = pj . 
Therefore, Eq. ([6]) is a concave quantity, so that ineq. (jH]) 
is indeed convex. Since, as shown above, it is satisfied for 
all biseparable pure states, this implies the same also for 
mixed states. 

Ineqs. HI and (|TTl are valid for any choice of a com- 
pletely separable pure state- vector |$), but the potential 
to detect the genuine multipartite character of a given 
entangled state will depend on a suitable choice of |$). 
For a state with rather weak genuine multipartite entan- 
glement an optimization, i.e. a search for the product- 
vector 1$) that maximizes the violation of the respective 
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FIG. 1: Here the detection quality of the bilinear inequali- 
ties dU, (mi and (fill)) is shown for the state p = + 
apGHZ + /3pw (a) and the tripartite qutrit state (with sub- 
systems labeled ABC) p = ^—^^1 + aptiaep + PpgGHZ (b). 
Area [nj contains genuine multipartite entanglement detected 
by (|lT]|. Area IIIll contains genuine multipartite entanglement 
detected by (jlTI]). Area|T|(a) is not biseparable w.r.t. any 
bipartition, since it violates inequality (jl]) for all partitions, 
area|T]i3|j4C is not biseparable w.r.t. B\AC, since it violates 
inequality ((I)) for this partition (the result for AB\C is equiv- 
alent). Area|llA|_BC contains states that violate inequality 
^ for bipartition A\BC. The area labeled PPT constitutes 
all states not detected by the Peres-Horodecki criterion '5]. 



inequality might be necessary. Optimizations that are 
quite common in the theory of entangled states, and gen- 
erally pose a difficult problem. But, finding an optimal 
state vector is significantly easier than e.g. convex roof 
construction or optimization of entanglement wit- 
nesses since there are no constraints to be satisfied 
and efficient algorithms are available 14|. However, what 



is even more important, the optimization space grows 
only linearly with the number of subsystems as opposed 
to the typically exponential scaling of such problems. Be- 
sides such an optimization, one can combine different 
choices of states |<I>) to tailor criteria that are suited par- 
ticularly well for a specific class of states, as we demon- 
strate here with the exemplary choice of = \si) ®\sj) 
with I Si) = \x . . . xyx ... a;) in terms of two single-particle 
states \x) and jy), and \y) is chosen exactly for the i-ih 
entry of \si). Taking linear combinations of ineq. for 
these choices we arrive at 



^ ^($,,|p®2n|<i>,^.) - (n - 2) ^ ^{^^Mp'^^m^,) < 

(III) 

where Vi — IlAi^'i-Bi is defined analogously to the above. 
However, in contrast to the above, not all bipartitions 
are taken into account, but Ai is the duplicated Hilbert 
space of the i-th subsystem and Bi the rest. Exactly as in 
ineq. (jlTl also the left-hand-side in ineq. (jIII[) is convex, so 
that the inequality is proven for biseparable mixed states, 
since it is proven for biseparable pure states in the end 
of the paper. 

In particular with growing system size the ability to 



assess separability criteria efficiently is getting more and 
more important since quantum state tomography (QST) 
scales so unfavorably with the number of subsystems. 
Being bilinear expectation values, the present criteria 
can very efficiently be measured experimentally on iden- 
tically prepared quantum states as it has been done in 
[l5| . But also with measurements performed on individ- 
ually prepared quantum states the present criteria can be 
experimentally assessed with significantly fewer observ- 
ables than required for QST: Eq. pl| is given in terms 
the square root of the number of observables needed for 
QST and Eq. (IIIII) scales as 2n(n — 1), i.e. polynomially 
with the number n of subsystems. 

An indispensable prerequisite for any practical crite- 
rion is its robustness against experimental imperfections. 
Therefore, let us discuss the capacity of our present cri- 
teria on a few exemplary quantum states, where we test 
this robustness. 

Example 1. First consider the three qubit state p = 
"^~g~^ l -I- apGHZ + P pw where pcHZ = 1(1000) + 
|111))((000| + (111|) and pw = jdOOl) + |010) + 
100))((001|-|-(010|-|-(100|). It is a mixture of the GHZ- 
state and the W^-state dampened by isotropic noise (see 
Ref.flil for further details). In Fig. [T^ the detection pa- 
rameter spaces of the inequalities (HI, (|n| and (jIII[) are 
illustrated. In the case of genuine multipartite entangle- 
ment detection for qubits, these criteria work as well as 
the best known method so far. For example in Ref. [13] 
the above state for (a = and f3 = 1 — p) was found to 
be genuinely multipartite entangled by means of entan- 
glement witnesses up to a threshold ofo < 8/19. This 
bound was then improved to p < 8/17 [18| . which is also 
our result. In fact for this special case our criteria co- 
incide. For qudits, our criteria are the first detection 
criteria known so far. 

Example 2. Consider the three qutrit state p — ^—^^^\~\- 
+ PpgGHZ where pgCHZ = sdOOO) -I- |111) + 
|222))((000|-H(lll|-|-(222|)and pusep = i(|0)(0|«)(|00) + 
|11) + |22))((00| + (11| + (22|). It is a mixture between 
a generalized GH Z-staic for qutrits and a biseparable 
qutrit state dampened by isotropic noise. In Fig. [TJd 
the detection parameter spaces of the violation of the 
inequalities ^ and (jlll are illustrated, 
g Example 3. Now consider the following four qudit state: 



PS 



d4 



- 1 + ^ E f^sGHZl + ^ E PIgHZ2 (8) 



where ffgGHZx ■= \gGHZx{i)){gGHZx{i)\ with 
\gGHZx{i)) := J2k ^|A:)|fc®a;)|A:©i)|fc©i©a;). Where 
® is the addition modulo d. For d — 2 and a = /3 this is 
the bound entangled Smolin state (see Ref. [l^) damp- 
ened by isotropic noise. Also in this case our criteria 
work well. Ineq. ([3]) shows that all states in the region 
l-{d'^ + l)a-P < and l-a-{d'^ + l)l3 < are not sepa- 
rable with respect to any bipartition. Moreover ineq. pi)) 
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shows that the state only becomes biseparable outside the 
region not detected by ineq. ([3]) for a = /3 > l/{d? + 2), 
i.e. all entangled states in this region except for the line 
a = (3 are definitely multipartite entangled. That is, our 
criteria detect all states that are detected by the Peres- 
Horodecki criterion but, - exceeding the scope of the 
latter - characterize states to be genuinely multipartite 
entangled. 

Example 4 Consider the generalized GHZ state {ipdn) = 
^ J2iZo K)*^" with additional isotropic (white) noise: 



P = Pli'dn) {ijMn \ + (1 



(9) 



With inequality (|TI]) we can show analytically that these 
states are genuinely multipartite entangled for p > 
d"-i+3 ' which shows that even in high dimensional sys- 
tems with many constituents these criteria work very 
well. 

In conclusion ineq. (jIII[) is only one specific of many 
possible criteria derived from ineq. pB and the versatility 
of our approach allows to tailor many criteria suited for 
specific classes of entangled states. Given the efficient 
decomposition in physical observables, our criteria enable 
verification of mixed state entanglement with tools (20l . 
[2H originally only applicable to pure states. 
Appendix: Finally, let us prove ineqs. (j3]) and For 
the former, we have to show 

for any separable mixed state Qs = X^i I ^ 

\Xi){Xi\ and any completely separable state- vector |<i>) — 
® This amounts to showing 



X*X - 



with [X]pi...p„qi. 

[^AX]pi...Pmqi---q7i 
Since (n^/sX)* (iIa/bX 



UbX 

nr=i 



fies to 



UaX - UbX 



- - (UbX) ' (uax) > , 
(10) 

I\T=l{(^■^\'Pp^)IlT=l{^^i\xq^), 

"*l'^P.+imod™^nili(ftlX9.), 

X*X, ineq. (HU]) simpli- 
> 0, which proves ineq.®. 



Similar to the proof presented in Ref. [iSj we 
only have to verify that ineq. (jllip is satisfied for 
pure biseparable states |^) due to its convexity (as 
shown for ineq. ^Qi). With the short hand notation 

X.,, - ^($,,|pfn|$,,) and y,, - ^($,, ^,1$,,), 

ineq. (|lll| reads J^i^j^i] - [n - '^)Y,ijVi]- We will 
have to distinguish between the cases in which both in- 
dices i and j correspond to different, or the same parts 
A and B in the bipartition with respect to which \^) 
(without loss of generality we assume i to correspond 
to A). The former contributions to ineq. (IIIII) we de- 
note as Bd = ^ij^B^^ij — [n — '2)yij), the latter as 



{n - 2)Y,^yit, so 
Bd is non-positive 
Sh 



Bs = T,^^,eAx^3 - {n - 2)yy) - 
that ineq. (jnll reads B^ + Bd < 0. 
since Xij < i/ij as shown for ineq. (|Tl|. Since the yij 
are non-negative, we obtain Bg < J2i^j£s — {n — 
2) J2^ y^^ = J2^^jesi^^] " ^^y«) < E,;#iGsfej - with 
Zi — [n ~ 2)/(ni — 1), where rii is the number of sub- 
systems in A, where Zi > 1 since A can comprise at 
maximum n — 1 subsystems. Now, we can symmetrize 
the last term in '^^^j^g{xij — ya), i.e. rewrite it as 
I]i/jes(^ii ~ 1/2(2/^ + y]])). Since yu = {si\Q-q,\si) 
(due to the relation Vi\^ii) — we can conclude 

Xij = |(si|£i*|si)| < ll2{yii + yjj)), such that Bg is non- 
negative, what finishes the proof of ineq. (|IIip . 
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